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Abstract 

We present a mean-field theory of carrier-induced ferromagnetism in Mna;IIIi_2,V diluted magnetic semiconductors 
with a special emphasis placed on the magnetic anisotropy. The valence band holes are described using the six 
band Kohn-Luttinger model. We find that the magnetic anisotropy is a complicated function of sample parameters 
such as hole density or strain. Results of our numerical simulations are in agreement with magnetic anisotropy 
measurements on samples with both compressive and tensile strains. 
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1. Introduction 

Experiments [1,2] in Mna;IIIi_2;V diluted mag- 
netic semiconductors (DMS) have demonstrated 
that these ferromagnets have remarkably square 
hysteresis loops with coercivities typically ~ 40Oe, 
and that the magnetic easy axis is dependent on 
epitaxial growth lattice-matching strains. The 
physical origin of the anisotropy energy in our 
model is spin-orbit coupling in the valence band. 
Our work is based on six band envelope function 
description of the valence band holes and a mean- 
field treatment of their exchange interactions 
with Mn"*"*" ions. Even in the mean-field theory, 
we find that the magnetic anisotropy physics of 



these materials is rich and we predict easy axis 
reorientations as a function of sample parameters 
including hole density or epitaxial growth lattice- 
matching strains. Similar conclusions have been 
presented in a closely related, independent study 
[3]. A formal theory of magnetic anisotropy in 
cubic semiconductor DMS is derived in Section 2. 
Numerical results for typical experimental samples 
are presented in Section 3. 
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2. Formal theory 

In zero external fields the Hamiltonian for the 
valence band electrons interacting with localized d 
electrons on the Mn++ ions can be written as [4] 
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where i labels a valence band hole and / labels a 
magnetic ion, S/ is a localized spin, Sj is a hole spin, 
and Hl is a six-band envelope-function Hamilto- 
nian [5] for the valence bands. In IIIi-a^Mnj^V semi- 
conductors, the four j = 3/2 bands are separated 
by a spin-orbit splitting Ago from the two j = 1/2 
bands. In the relevant range of hole and Mn"'""'" 
densities, no more than four bands are ever occu- 
pied. Nevertheless, mixing between j = 3/2 and 
j = 1/2 bands does occur, and it can alter the bal- 
ance of delicate cancellations which often controls 
the net anisotropy energy. The exchange interac- 
tion between valence band holes and localized mo- 
ments is believed to be antiferromagnetic [4], i.e. 
Jpd > 0. For GaAs, experimental estimates [6,7] of 
Jpd fall between 0.04eVnm^ and 0.15eVnm^, with 
more recent work suggesting a value toward the 
lower end of this range. 

In the mean-field theory [8,9] the valence band 
hole experience an effective field which depends on 
the average magnitude, M, and orientation, M, of 
the localized moments: 

hMF(M) = JpdNMnMM = hM (2) 

where Nmu = Nj/Vis the number of magnetic im- 
purities per unit volume. The magnetic anisotropy 
in the absence of strain is well described by a cu- 
bic harmonic expansion truncated at sixth order, 
an approximation commonly used in the litera- 
ture [10] on magnetic materials. The correspond- 
ing cubic harmonic expansion for total energy of 
a system of non-interacting fermions with single- 
particle Hamiltonian Hl — hM ■ s is 

Eb{M) = £6((100)) + Kl^MlMl + M^Ml 

+ MlMl) + K^°- M^M^M^. (3) 

The cubic anisotropy coefficients are related to to- 
tal energies for M along the high symmetry crystal 
directions by following expressions: 

- V 

_ 27^b(lll) - 36E,{nO) + 9^6(100) 
^2 - y -W 



The MBE growth techniques produce IIIi_a;Mna; V 
films whose lattices are locked to those of their 
substrates. X-ray diffraction studies [11] have es- 
tablished that the resulting strains are not relaxed 
by dislocations or other defects, even for thick 
films. Strains in the IIIi_2;Mn2;V film break the 
cubic symmetry assumed in Eq. (3). However, the 
influence of MBE growth lattice-matching strains 
on the hole bands of cubic semiconductors is well 
understood [12] and we can use the same formal 
mean-field theory to account for strain effects on 
magnetic anisotropy. 



3. Numerical results 

We turn now to a series of illustrative cal- 
culations intended to closely model the ground 

state of Ga.95Mn.05As. For this Mn density and 
the smaller values of Jpd favored by recent es- 
timates, h = Jpd^MnJ ~ O.OlRy at zero tem- 
perature. This value of h is not so much smaller 
than the spin-orbit splitting parameter in GaAs 
(Ago = 0.025 Ry), so that accurate calculations 
require a six band model. Even with x fixed, our 
calculations show that the magnetic anisotropy of 
Ga.95Mn.Q5As ferromagnets is strongly dependent 
on both hole density and strain. The hole density 
can be varied by changing growth conditions or by 
adding other dopants to the material, and strain 
in a Ga.95Mn.05As film can be altered by changing 
substrates. The cubic anisotropy coefficients (in 
units of energy per volume) for strain- free material 
are plotted as a function of hole density in Fig. 1. 
The easy axis is nearly always determined by the 
leading cubic anisotropy coefficient Kf^, except 
near values of p where this coefficient vanishes. As 
a consequence the easy-axis in strain free samples 
is almost always either along one of the cube edge 
directions {Kf^ > 0) , or along one of the cube di- 
agonal directions (JsTf" < 0). Transitions in which 
the easy axis moves between these two directions 
occur twice over the range of hole densities stud- 
ied. (Similar transitions occur as a fimction of h, 
and therefore temperature, for fixed hole density.) 
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Fig. 1. Cubic magnetic anisotropy coefficients as a function 
of hole density 

Near the hole density 0.01 nm^'^, both anisotropy 
coefficients vanish and a fine-tuned isotropy is 
achieved. The slopes of the anisotropy coefficient 
curves vary as the number of occupied bands in- 
creases from 1 to 4 with increasing hole density. 
This behavior is clearly seen from the correlation 
between oscillations of the anisotropy coefhcients 
and onsets of higher band occupations. 

Six-band model Fermi surfaces are illustrated in 
Figs. 2 and 3 by plotting their intersections with 
the kz = plane at p = O.lnm"'^ for the cases of 
(100) and (110) ordered moment orientations. The 
dependence of quasiparticle band structure on 
ordered moment orientation, apparent in compar- 
ing these figures, should lead to large anisotropic 
magnetoresistance effects in llli-^MnxV ferro- 
magnets. We also note that in the case of cube 
edge orientations, the Fermi surfaces of different 
bands intersect. This property could have impor- 
tant implications for the decay of long- wavelength 
collective modes. 

In Fig. 4 we present mean-field theory predic- 
tions for the strain-dependence of the anisotropy 
energy a,t h — O.OlRy and hole density p = 
0.35nm~^. According to our calculations, the easy 
axes in the absence of strain are along the cube 
edges in this case. This calculation is thus for a 
hole density approximately three times smaller 
than the Mn density, as indicated by recent exper- 
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Fig. 2. Six-band model Fermi surface intersections with 
the fcz = plane for p = O.lnm"'^ and h = O.OlRy. This 
figure is for magnetization orientation is along the (100) 
direction. 
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Fig. 3. Six-band model Fermi surface intersections with the 
fcz = plane for the parameters of figure 2 and magneti- 
zation orientation along the (110) direction. 
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strains, the sign of the anisotropy changes empha- 
sizing the subtlety of these effects and the latitude 
which exists for strain-engineering of magnetic 
properties. 
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Fig. 4. Energy differences among (001), (100), (110), and 
(111) magnetization orientations vs. in-plane strain eo at 
h = 0.01 Ry and p = 3.5 nm"''. For compressive strains 
(eo < 0) the systems has an easy magnetic plane perpen- 
dicular to the growth direction. For tensile strains (eo > 0) 
the anisotropy is easy-axis with the preferred magnetiza- 
tion orientation along the growth direction. The anisotropy 
changes sign at large tensile strain. 

iments. The relevant value of the in-planc strain 
produced by the substrate-film lattice mismatch, 



eo 



as - af 
a/ 



(5) 



depends on the substrate on which the epitaxial 
Ga.95Mn.05As film is grown. The most important 
conclusion from Fig. 4 is that strains as small as 
1% are sufficient to completely alter the magnetic 
anisotropy energy landscape. For example for 
(Ga,Mn)As on GaAs, eo = -.0028 at a; = 0.05, the 
anisotropy has a relatively strong uniaxial contri- 
bution even for this relatively modest compressive 
strain, which favors in-plane moment orienta- 
tions, in agreement with experiment. A relatively 
small (^1 kJ m~'^) residual plane-anisotropy re- 
mains which favors (110) over (100). For x — 0.05 
(Ga,Mn)As on a x 0.15 (ln,Ga)As buffer the 
strain is tensile, eo — 0.0077, and we predict a 
substantial uniaxial contribution to the anisotropy 
energy which favors growth direction orientations, 
again in agreement with experiment. For the ten- 
sile case, the anisotropy energy changes more dra- 
matically than for compressive strains due to the 
depopulation of higher subbands. At large tensile 
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